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Abstract. In this paper, we introduce polynomial time algorithms that generate random fc- 
noncrossing partitions and 2-regular, fc-noncrossing partitions with uniform probability. A fc- 
noncrossing partition does not contain any k mutually crossing arcs in its canonical representation 
and is 2-regular if the latter does not contain arcs of the form + 1). Using a bijection of 
Chen et al. [2j [4] , we interpret fc-noncrossing partitions and 2-regular, fc-noncrossing partitions as 
restricted generalized vacillating tableaux. Furthermore, we interpret the tableaux as sampling 
paths of a Markov-processes over shapes and derive their transition probabilities. 



1. Introduction 

Recently Chen et. al. [3] studied a computational paradigm for fc-noncrossing RNA pscudoknot 
structures. The authors regarded fc-noncrossing matchings as oscillating tableaux [2] and inter- 
preted the latter as stochastic processes over Young tableaux of less than fc rows. According 
to which, they generated fc-noncrossing RNA pseudoknot structures [5] with uniform probability. 
Since the generating function of oscillating lattice walks in Z fc_1 , that remain in the interior of 
the dominant Weyl chamber is given as a determinant of Bessel functions [5], it is D- finite. As a 
result, the transition probabilities of the Markov- process can be derived in polynomial time as a 
pre-processing step and each fc-noncrossing RNA pseudoknot structure can be generated in linear 
time. 



The analogue of the above result is for fc-noncrossing partitions much more involved. Only for k = 3, 
Bousquet-Melou and Xin [T] compute the generating function and conjecture that fc-noncrossing 
partitions are not P-recursive for k > 4. Since the enumerative result are unavailable for fc > 4, it 
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is important to construct an algorithm which can uniformly generate fc-noncrossing partitions in 
order to obtain statistic results. 

In this paper, we show that fc-noncrossing partitions and 2-regular, fc-noncrossing partitions can be 
uniformly generated in linear time. We remark that there does not exist a general framework for 
the uniform generation of elements of a non-inductive combinatorial class. In this context relevant 
work has been done by Wilf [9] as well as [7] . 

The paper is organized as follows: in Section [2j we present all basic facts on fc-noncrossing par- 
titions, vacillating tableaux and lattice walks. In Section [3j we compute the number of lattice 
walks ending at arbitrary v £ Qk-i from the construction and prove Theorem 13. 2\ the uniform 
generation result for fc-noncrossing partitions. In Section [5J we show how to uniformly generate 
2-rcgular, fc-noncrossing partitions. To this end, we establish a bijection between lattice walks as- 
sociated to partitions and braids [4] . We will show that walks associated with 2-regular partitions 
correspond to walks associated to loop-free braids. The latter can easily be dealt with via the 
inclusion-exclusion principle. In Theorem 14.41 wc present an algorithm that uniformly generates 
2-regular, 3-noncrossing partitions in linear time. Finally, in Section [5J we show the particular 
algorithms in case of fc = 4. 

Wc remark that the generation of random 2-rcgular, fc-noncrossing partitions has important ap- 
plications in computational biology. The latter allow to represent so called base triples in three 
dimensional RNA structures, see [8] for details of the framework. The uniform generation then 
facilitates energy-based, ab initio folding algorithms along the lines of [3] . This paper is accompa- 
nied by supplemental materials containing MAPLE implementations of the algorithms, available 
at www.combinatorics.cn/cbpc/par.html. 



2. Some basic facts 

A set partition P of [N] = {1,2,..., N} is a collection of nonempty and mutually disjoint subsets 
of [N], called blocks, whose union is [N]. A fc-noncrossing partition is called m-regular, m > 1, 
if for any two distinct elements x, y in the same block, we have \x — y\ > m. A partial matching 
and a matching is a particular type of partition having block size at most two and exactly two, 
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respectively. Their standard representation is a unique graph on the vertex set [N] whose edge set 
consists of arcs connecting the elements of each block in numerical order, see Fig[TJ 

maybe new similar graph 



FIGURE 1. Standard representation and fc-crossings: we display the partition {2}, {8}, 
{11}, {1,5}, {6,12}, {9,15}, {13,14}, {3,7,10} of [15]. Elements within blocks are 
connected in numerical order. The arcs {(1, 5), (3, 7)}, {(3, 7), (6, 12)} {(7, 10), (9, 15)} 
and {(6, 12), (9, 15)} are all 2-crossings. 



Given a (set) partition P, a k-crossing is a set of fc edges (*2,j2), ■ ■ • , (*fe, jk)} such that 

i\ < ii < . . . < ik < ji < ji < ■ . ■ < jki sec FigfT] A k-noncrossing partition is a partition without 
any fc-crossings. We denote the sets of fc-noncrossing partitions and m-regular, fc-noncrossing 
partitions by CPfc(iV) and J , fc,m(^V), respectively. For instance, the set of 2-regular, 3-noncrossing 
partitions are denoted by J > 3,2(A f ). 

A (generalized) vacillating tableau [4] V£ N of shape A and length 2N is a sequence A , A 1 , ... , A 2JV 
of shapes such that (1) A = 0, \ 2N = A and (2) for 1 < i < N, A 2 * -1 , A 2 * are derived from A 2 *~ 2 
by elementary moves (EM) defined as follows: (0,0): do nothing twice; (— D,0): first remove a 
square then do nothing; (0, +□): first do nothing then add a square; (±D,±D): add/remove a 
square at the odd and even steps, respectively. We use the following notation: if Ai+i is obtained 
from A; by adding, removing a square from the j-th row, or doing nothing we write A,+i \ A; = +□■/, 
Aj+i \ Aj = —Dj or Aj-|_i \ A, = 0, respectively, see Figf2] 

A braid over [N] can be represented via introducing loops and drawing arcs and 
with i < j < £ as crossing, see Fig(3](Bl). A k-noncrossing braid is a braid without any /c-crossings. 
We denote the set of fc-noncrossing braids over [TV] with and without isolated points by 53^ (N) 
and 23£(AT), respectively Chen et al. [2] have shown that each fc-noncrossing partition corresponds 
uniquely to a vacillating tableau of empty shape, having at most (fc — 1) rows, obtained via the 
EMs {(-111,0), (0,+D), (0, 0), (-□,+□)}, see Figf3](A2). In [4], Chen et al. proceed by proving 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 
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(A2) :<f> +□: +n i -u'\t> :-□ +□!;*+□.. +□ -□.:'+□, 
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-□ !:-□ ! 



(C1) 00 □□□□□□ gang nng a a □ □□00 

(C2) !0 +oj j+o -□. i : i +□ -□ M +□ :!-□ 1: 1 -: : +i 1 -d^ jl-o 0; 

FIGURE 2. Vacillating tableaux and elementary moves: (Al) shows a general vacillating 
tableaux with EM (A2). (Bl) displays another vacillating tableaux and its set of EM, 
{(-□,0), (0,+D), (0,0), (-□,+□)} shown in (B2). In (CI) we present a vacillating 
tableaux with the EM {(-□, 0), (0, +□), (0, 0), (+□, -□)} displayed in (C2). 



that vacillating tableaux of empty shape, having at most (k — l)-rows which are obtained by the 
EMs {(— □, 0), (0, +□), (0, 0), (+□, — □)} correspond uniquely to fc-noncrossing braids. 



1 2 3 4 5 6 7 8 9 10 11 12131415 1 2 3 4 5 6 7 8 9 10 11121314 

(A2) 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 1920 212223 24 25 26 27 28 29 30 

00a □ □□ g g g □ □ □ g □ mmmmg: g g g g □ □ □□□□□00 

(B2) 

01 23 456789 10 11 12 13 14 15 16 17 18 1920 212223 24 25 26 27 28 

00nnnngnnnngnnmmmgng g g □ □ □ □ m 000 

FIGURE 3. Vacillating tableaux, partitions and braids: in (Al) we show a 3-noncrossing 
partition and in (A2) the associated vacillating tableau. (Bl) shows a 3-noncrossing braid 
and (B2) its vacillating tableau. 



In the following, we consider A:-noncrossing partition or a 2-rcgular, fc-noncrossing partition of TV 
vertices, where N > 1. To this end, we introduce two Zfc_i domains Qk-i = {v = (yi, ^2, ■ ■ ■ , ^fc-i) € 
Zfc_i I vi,V2, ■ • ■ , Vk-i > 0} and Wk-i = {v = (vx, ^2, ■ ■ ■ , Vk-x) G %>k-i | ^1 > ^2 > • ■ • > Vk-i > 
0}. Set Dk-i £ {Qk-i, Wk-i}, e-i be the unit of i-th axis, where i £ {1, 2, . . . , k — 1} and eo = 
be the origin. A Tr) k _ 1 -walk is a lattice walk in D^-i starting at e = {k — 2, , k — 1, . . . , 0) having 
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steps dbej, such that even steps are +ei and odd steps are either — ej, where i € {0, 1, — 1}. 
Analogously, a r B>r> kl -walk is a lattice walk in D^-i starting at e whose even steps are either — ei 
and whose odd steps arc +e,, where i E {0, 1, . . . , k — 1}. By abuse of language, we will omit the 
subscript D^-i- 

Interpreting the number of squares in the rows of the shapes as coordinates of lattice points, we 
immediately obtain 

Theorem 2.1. @] 

(1) The number of k-noncrossing partitions over [N] equals the number of fwu-i -walks from e to 
itself of length 2N. 

(2) The number of k-noncrossing braids over [N] equals the number oj "3w s ,_ 1 -wato from e to itself 
of length 2N. 



Recall 5 = (k — 2, , k — 1, . . . , 0), Bousquet-Melou and Xin p] according to the reflection principle, 
arrive at the following proposition: 

Proposition 2.2. For any ending points (i £ Wk—i> set uj^i and a ke denote the number ofWk—i- 
walks and Qk-i going from 5 to fi of length I, respectively. Then we have 

(2.1) Y, sg^K?' 

where is the permutation group of [k — 1], sgn(7r) is the sign of it and 

(2-2) 7r(/x) = 7r(^i,/! 2 , . . . ,/ife-l) = (Ai T (l),M7r(2), • ■ • ,Mir(fc-l))- 



3. RANDOM fc-NONCROSSING PARTITIONS 



In this section, we establish an algorithm to uniformly generate a fc-noncrossing partition of [N] 
employing a Markov-process as an interpretation of a vacillating tableau. 

Lemma 3.1. (a) Suppose 1 < £ < N — 1 and v = (i>\, ■ ■ ■ , Vk—\) G Qk-i, then set a v k 2i+1 = f% 
accordingly we have 



(3.1) 



f v h>l for s = 21 + 1, 

EtoXP for S - 2^ + 2. 
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where f v kl = for £\ i/ 4 > £ + (fc ~ 1} 2 (fc ~ 2) and /or £ = 0, 

fk,0 

(b) ffcg satisfies the recursion 

(3.2) /M=EECr ej 



1 /or ^ = A; 
otherwise. 



k-l k-l 



i=0 i=0 



Proof. Assertion (a) is obvious by construction. To prove assertion (b). indeed by definition of 
fk e> we nave a k 2 e+i = fk v Furthermore, since each lattice walk ending at v of length 21 + 1 can 
be obtained from a lattice walk ending at (v — G\ + ej) G Wk-i of length 2£ — 1 via first +ej then 
— ej. I.e. we obtain cq. 13.21 and the proof of the lemma is complete. □ 



In the following, we consider a vacillating tableau as the sampling path of a Markov-process, whose 
transition probabilities can be calculated via the terms ui% £■ 

Theorem 3.2. A random k-noncrossing partition can be generated, in 0(N k ) pre-processing time 
and 0(N k ) space complexity, with uniform probability in linear time. Each k-noncrossing partition 
is generated with O(N) space and time complexity, see Algorithm 1. 



Proof. The main idea is to interpret tableaux of fc-noncrossing partitions as sampling paths of a 
stochastic process, see FigH]for the case k = 3. We label the (i + l)-th shape, A^ +1 by a = A l+1 \ A 4 
and a G — Qy, 0}jZ 1 -1 , specifying the particular transition from A' to A i+1 . 

Let Vfc(A„ +1 , 2N — (i + 1)) denote the number of vacillating tableaux of length (i + 1) such that 
A i+1 \ A' = a. We remark here that if A i+1 has (vi — k — i + 3) boxes in the i-th row, then 
V fc (AJ +1 ,2A^- (i + 1)) = ^k,2N-(i+i)' wn ere v = (ui, u 2 , . . . , ffc-i). Let (X 1 )^ be given as follows: 

• X° = X 2N = and X 1 is a shape having at most k — 1 rows; 

• for 1 < i < N - 1, we have X 2i+1 \ X 2l e {0, -D^Z^ 1 and X 2l+2 \ X 2t+1 G {0, +n j }^~ 1 . 

• for 1 < i < 2N — 1, we have 



(3.3) P(X l+1 = A l+1 | X 1 = A 4 ) = V(A " +1 ' 2N 1 l) 



V k (A*,2JV-i) 
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Algorithm 1 Uniform generation of fc-noncrossing partitions 
1: 1=1 

2: Tableaux (Initialize the sequence of shapes, {A l }-f 

3: A = 0, X 2N = 

4: while i < 2N do 

5: if i is even then 

6: X[0]«- V k (Ag- x , 2JV - (i + 1)) 

7: for 1 < j < k — 1 do 

8: X[j]^-V k (A^,2JV-(i + l)) 

9: end for 

10: end if 

11: if i is odd then 

12: X[0]^V k (A^ +1 , 2JV-(i + l)) 

13: for 1 < j < k — 1 do 

i4: x[i]^v k (A!+4,2J\r-(i + i)) 

15: end for 
16: end if 

17: ^m^(E^i _1 XLj])+X[0] 

18: Shape <— Random(sum) (Random generates the random shape ^p./ or _p.) with probability 
X[j]/ sum or A^ 1 with probability X[0]/sum) 

19: i <— i + 1 

20: Insert Shape into Tableau (the sequence of shapes). 
21: end while 

22: Map( Tableau) (maps Tableau into its corresponding fc-noncrossing partition) 



In view of eq. (|3 . 3[) . we immediately observe 



2AT-1 



(3.4) [ P(X 1+1 = A l+1 | X 1 = A 4 ) 



i=0 



V k (A^=0,O)_ 1 
V k (A° = 0,iV) V k (0,2iV)' 



Consequently, the process (X z )f^ generates random fc-noncrossing partitions with uniform prob- 
ability in O(N) time and space. 

As for the derivation of the transition probabilities, suppose we are given a shape X h having exactly 
(i>i — k — i + 3) boxes in the i-th row. where i G {1, 2,...,k— 1}. According to Lemma |3.1[ set 
v = (i>i, V2, ■ ■ ■ , Vk-\)i then we reduce the problem obtain to^ to the calculation of several coeffi- 
cients /j? for some fixed (3 £ Wk-i related to v and 1 < £ < N — 1. 
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0123 4 5 678 

m -* □ 

cfi — — □ —►□—»- □ ^ — 

^ □ — □ —0 ► 

V g -» □ — □ — —0 

FIGURE 4. Sampling paths: we display all tableaux sequences where A 1 = 0, A 2 = 
A 3 = □. We display the three possible A 4 -shapes (blue), induced by , +O2 and 
all possible continuations (grey). 



Furthermore, according to Lemma [3TT1 (b). the coefficients fg for all (3 G Wk-i and 1 < £ < N — 1, 
can be computed in 0(N k ) time and 0(N k ) space complexity via fc-nested For-loops. 
Given the coefficients ff for all (3,1, we can derive the transition probabilities in 0(1) time. Ac- 
cordingly, we obtain all transition probabilities V|<(A^,27V — i)) in 0(N k ) time and 0(N k ) space 
complexity. □ 



4. RANDOM 2-REGULAR, fc-NONCROSSING PARTITIONS 



In this section, we generate random 2-regular, fc-noncrossing partitions employing a bijection be- 
tween the set of 2-regular, fc-noncrossing partitions over [N] and the set of fc-noncrossing braids 
without loops over [N — 1]. 

Lemma 4.1. Set k G N, k > 3. Suppose (X 1 ) 2 ^ 1 in which X° = is a sequence of shapes such 
thatX 2 ^\X 2 ^ 1 G {{+a h }', l = k -\0} andX^-^X 2 !- 2 G {{-□ /l }^Z^ 1 , 0}. Then (X 1 ) 2 ^ 1 induces 
a unique sequence of shapes (/^)f£ with the following properties 

(4.1) \ ^ G {{+U h } h h Z\-\0} and ^+ 2 \ G {{-n h }£*-\ 0}, 

(4.2) ^ = X 2 > +1 

(4.3) fi 2j+1 ^ X 2j+2 =► m 2j+1 S {\ 2j+1 ,X 2 i+ 3 }. 
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partition 




11 12 13 14 15 



Vacillating 1 2 3 4 5 6 7 
tableaux 



9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 
sequence <& <fi □ □□ □ B B ^ B n ^ rn rn |TI |T| Q j~j j~j j~j fl pi 11 □ □ □ ^6 



<p □□ □□ g BBB B D B nmn=||=|=|n=| =| nnnni=l □ 



Elementary 

move 

sequence 



00 



+P«6i 



0; 



,60 



00 



* * 0-n ! 



Elementary 

move 

sequence 



0+D 



00 



0+D, 







□ +□ -C -C 



00 



00 



-□ 0: 







□ 



tin g □ □ □ □ gggnagnnmmmrpgggnnma □ n i> <f> 



12 3 4 



9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 



3-noncrossing 
braid -] 




10 11 12 13 14 



FIGURE 5. Illustration of Lemma l4.ll We show how to map the ^Py^-walk induced by 
the 3-noncrossing partition P, of Fig|3] (A) into the 25 w 2 -walk. The latter corresponds 
to the 3-noncrossing braid of Fig(3] (B). Here each vertex, i, is aligned with the triple of 
shapes (A 2 ' -2 , A 2 ' -1 , A 2 ') in the corresponding tableaux. 



Proof. Since A 1 = 0, (A- 7 ) 2 ^ 1 corresponds to a sequence of pairs ((xi, yi)) l i=l given by 

X 2i \ X 2i-1 and y , = A 2*+l \ X 2i such that 

(4.4) Vl<z<^; faM) e {(0,0), (+D h , 0), (0, -□/>), (+□&,-□,)}. 

Let </? be given by 



(4-5) v((a*,l«)) 



(a:i,2/i) for (a;,-, y*) = (+□?,, -Bf) 
{ViiXi) otherwise, 



and set (<p( Xi , Vl )) l l=1 = {a^b.^Li- Note that (a,, b t ) e {(-O h , 0), (0, +D h ), (0, 0), -Dj)}, 
where 1 < < — 1, see Fig. [SJ Let (//)f =1 be the sequence of shapes induced by (aj,6j)f_ 1 
according to a, = /i 2l \/i 2l_1 and 6j = /i 2 * -1 \^i 2l_2 initialized with /io = 0- Nowcq. (|4.ip is implied 
by eq. (|4.4p and eq. (|4. 2|) follows by construction. Suppose ^ 2j+1 ^ X 2 ^ +2 for some < j < £ — 1. 
By definition of ip, only pairs containing "0" in at least one coordinate arc transposed from which 
we can conclude /j 2j+1 = /i 2 -? or ^i 2j+1 = /i 2:,+2 , whence eq. ()4.3|) . I.e. we have the following 
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situation 



A 2 ' 



+1 



,2j + l 



,2j+2 



X 2j+3 



A 2i+1 



A 2j+3 




or 




/i 2 ^ /i 2 ^ 1 — A t 2 J+ 2 , 



and the lemma follows. 



□ 



Lemma FOl establishes a bijection between CPvi4_i-walks of length 2£ + l and B^^-walks of length 
21, where W k -i = {(01,02, . . . , a k -i) | Oi > a% > ■ ■ ■ > a fc _i}. 

Corollary 4.2. Let 1Pk(N) and r B k (N — 1) denote the set of k-noncrossing partitions on [N] and 
k-noncrossing braids on [N — 1] . Then 
(a) we have a bijection 



(4.6) 



1?: V k {N) — ►2 fc (i\T-l). 



(b) induces by restriction a bijection between 2-regular, k-noncrossing partitions on [N] and 
k-noncrossing braids without loops on [N — 1]. 



Proof. Assertion (a) follows from Lemma 14.11 since a partition is completely determined by its 
induced CPw^-walk of shape A 2 "^ 1 = 0. (b) follows immediately from the fact that, according to 
the definition of ip given in Lemma |4~T1 any pair of consecutive EMs (0, +Di), (— Di, 0) induces an 
EM (+Di, —di). Therefore, $ maps 2-regular, fc-noncrossing partitions into fc-noncrossing braids 
without loops. We illustrate Corollary |4.2f b) in Fig[51 □ 

Lemma 4.3. The number of 23^ -walks ending at v € W k -i of length 21, where 1 < i < N is 
given by 

(4-7) ^ = E(- 1 ) h Q w ^- h)+1 - 

Furthermore, we obtain the recurrence of o^'m-i 9^ ven by 
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.? SI 



+□,0; :00 



0+cr .00 



!00 



|00 



0-C 



0i 



■□ 0: 
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0-P.i 



□ 
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tableaux 
sequence 
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braid -] 




12 13 14 



FIGURE 6. Mapping 2-regular, 3-noncrossing partitions into loop-free braids: illustra- 
tion of Corollary I4.2f b). 



Proof. According to Lemma B~T1 any 3 3 vi/ fc _ 1 -walk of length 2£+l corresponds to a unique r Bw k _ 1 - 
walk of length 21, whence c£ 2 ^ = ^ 2 i+i- Let A2i(h) denote the set of 23 w k -± -walks of length 21 
in which there exist at least h pairs of shapes (fi2q-i, l^2q) induced by the EM (+ei, — ei). Since 
the removal of h EMs, (+ei, — ei), from a 25n/ fc _ 1 -walk, results in a 23w /i _ 1 -walk of length 2(1 — h), 
we derive A2t(h) = (f) ^hit—th- Using the inclusion-exclusion principle, we arrive at 

( 4 -9) a- = £(-!)* (fk, a( *-,o = E(- 1 ) h ff) w M(/- h)+1 - 



By construction, an odd step in a l &w h _ 1 -waXk is +ej, where j € {0, 1, . . . , k — 1}, whence 



k-l 

u:* _ V - " 1 - v— ej;* 
fc,2£-l ~~ 2 < a k.2e-2- 

3=0 



□ 
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Via Lemma [4.31 we have explicit knowledge about the numbers of 2 ^ fc _ 1 -walks, i.e. ay* for all 
v e Wk-i and 1 < s < 2N. Accordingly, we are now in position to generate 2-regular, k- 
noncrossing partitions with uniform probability via fc-noncrossing, loop- free braids. 

Theorem 4.4. A random 2-regular, k-noncrossing partition can be generated, in 0(N k+1 ) pre- 
processing time and 0(N k ) space complexity, with uniform probability in linear time. Each 2- 
regular, k-noncrossing partition is generated with O(N) space and time complexity, see Algorithm 
2, 3. 



Proof. We interpret 2^- -walks as sampling paths of a stochastic process. To this end, we again 
label the (i + l)-th shape, A 4 Q +1 by a = A i+1 \ A\ a £ {i+Oj}^' 1 , 0}, where the 

labeling specifies the transition from A' to A' +1 . In the following, we distinguish even and odd 
labeled shapes. 

For i = 2s, we let W£(\ 2 a s ,2N - 2s) denote the number of -walks such that a = A 2s \ A 25 " 1 . 

By construction, assume X 2s has i>i — k — i + 3 boxes in i-th row, respectively, then we have 



(4.10) 



WZ(\ 2 a s ,2N-2s) 



'k,2N-2s 



i.e. W^*(A 2s , 2N - 2s) is independent of a and we write W fe *(A 2s , 2N - 2s) = W£{\ 2s , 2N - 2s). 
For i = 2s + 1, let V£(\ 2 a s+1 , 2N - (2s + 1)) denote the number of 2^ _ t -walks of shape A 2s+1 of 
length 27V - (2s + 1) where A 2s+1 \ A 2s = a. Then we have setting u = 2N - 2s - 2 
(4.11) 



2s+l 



in which, 



e^oj w?(a^» + |(i + (-iy^(A 2 4+ 2 , U ) 



□i 

□ 



{WZ(\% s+2 ,u) 



for a — 
for a = 
for a = 0; 
for a = — □ 



j ■ 



J * 0, 1; 
3*0, 



(4.12) 



(3 = 



WZ(\ 2 a s+1 ,u+l) = Q, 
W* k {\ 2 a s +\u+l)^Q. 



We are now in position to specify the process (X 1 ) l ~^ N : 

• X° = X 2N = and X 1 is a shape with at most k — 1 rows 

• for l < i < n - 1, (x 2l+l \ x 2t , x 2l+2 \ x 2t ) e {(-□, 0), (0, +□), (0, 0), (+□, -□)} 

• there does not exist any subsequence (X 2 \ X 2i+1 , X 2i+2 ) such that (X 2l+1 \X 2 \ X 2i+2 \X 2i ) 

( + □!,-□!) 
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Algorithm 2 Tableaux: Uniform generation of 2-regular, fc-noncrossing partitions. 

1: Tableaux (Initialize the sequence of shapes to be a list {X l }iZo N ) 

2: A «- 0, X 2N «- 0, A 2Ar - J «- 0, i «- 1 

3: while i < 2N - 1 do 

4: if i is even then 

5: X[0]<- V£(A*+\2iV-(;+l)) 

6: for j = 1 to k — 1 do 

7: X[2j - 1]<- VJ(A^ ,2N-(i + 1)), X[2j]<- VJ(A^,2JV - (i + 1)) 

8: end for 

9: else 

10: X[0]<- \N*(X l +\2N - (i + 1)) 

11: for j = 1 to k — 1 do 

12: X[2j - 1]- WJ(A^ , 2iV - (i + 1)), X[2j]<- W k *(A!_+ n x . , 2iV - (i + 1)) 

13: end for 

14: FLAG(flagO, flagl, flag2, flag3) 

15: sum <- Etio" 1 ^^ 

16: end if 

17: Shape <— Random(sMm) (Random generates the random shape A^j .( or _Q .) with probability 

X[2j-l]/swm (or X[2j]/sum) or A^ 1-1 with probability X[0]/ sum) 

18: flagO <- 1, flagl <- 1, flag2 <- 1, flag3 <- 1 

19: if i is even and Shapc=A l f ^ i then 

20: flagl <— 

21: end if 

22: for 1 < s < k - 1 do 

23: if i is even and Shapc=A l+ [ -5 then 

24: flag2 «- 

25: end if 

26: end for 

27: for 2<s<fc-ldo 

28: if i is even and Shapc=A l f ^ then 

29: flag3 «- 

30: end if 

31: end for 

32: if % is even and Shapc=A^ hl then 

33: flagO <- 

34: end if 

35: Insert Shape into Tableaux (the sequence of shapes) 

36: i <— i + 1 

37: end while 

38: Map(Tableaux) (maps Tableaux into its corresponding 2-regular, fc-noncrossing partition) 
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Algorithm 3 FLAG: Distinguish the last odd step 

l: FLAG(ftagO, flagl, flag2, flag3) 

2: if flagO =0 then 

3: for j = 1 to k — 1 do 

4: X[2j]<- 

5: end for 

6: end if 

7: if flagl =0 then 

8: X[0]^ 0, X[2]<- 

9: for j = 1 to k — 1 do 

10: X[2j - 1]<- 

11: end for 

12: end if 

13: if flag2 =0 then 

14: for j = 1 to k — 1 do 

15: X[2j - 0, X[2j]<- 

16: end for 

17: end if 

18: if nag3 =0 then 

19: for j = 1 to k — 1 do 

20: X[2j - 1]<- 

21: end for 

22: end if 



• the transition probabilities are given as follows 

(1) for i = 21, we obtain 

P(X - A a |X - A ) - w , {xij2N _ i) • 

(2) for i = 21 + 1, we have 

p(x i + x X i + i lx i _ \i\ _ WZ(\ i+1 ,2N-i-l) 

nx A |X - AJ - v:(xi,2N-i) ■ 

By construction, 



2JV-1 



(4.13) [ P(X l+1 = A l+1 | X 1 = A*) = 



W fc *(A 27V =0,0) 



W fe *(A°=0,2iV) W fe *(0,2iV)' 
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whence (X l )f^ generates random 2-regular, fc-noncrossing partitions with uniform probability in 
O(N) time and space. 

As for the computation of the transition probabilities, according to Theorem 13 . 21 the terms co k e for 
v £ Wk-i and 1 < I < 2N + 1 can be calculated in 0(N k+1 ) time and 0(N k ) space. 
We claim that for fixed indices v = (y\, h>2, ■ ■ ■ , Vk-i) where v £ Wk-i and 1 < s < 2N, ov' can 
be computed in O(N) time. There are two cases: in case of s = 2£±, using a For-loop summing 
over the terms (— l)''^ 1 )^ 2(ti-h)+v we derive oy^- In case of s = 2i\ + 1, we calculate a v k ^'*, 
where % e {0, 1, . . . , k — 1} via cq. (|4.9| . Then oy 2^+1 follows in view of a v k 2^+1 — Sj=o a k 2I1* ■ 
Furthermore, using k nested For-loops for Vi, we derive a k J s for arbitrary v and s. Consequently, we 
compute for all v, s with 0{N k ) + 0(N k xJV) = 0(N k+1 ) time and 0(N k ) space complexity. 
Once the terms a,' for v and s are calculated, we can compute the transition probabilities in 
O(l) time. Therefore we obtain the transition probabilities W*(X l , 2N — i) in 0(N k+1 ) time and 
0(N k ) space complexity and the theorem follows. □ 



Remark We remark here it is feasible to generate a 2-regular, fc-noncrossing partition in 0(N k ) 
time complexity if we obtain the transition probabilities for 23^ fc -walks follow the similar rou- 
tine of -walks. The reason we show a different routine is we prefer to compute transition 
probabilities via the relation between two different combinatorial objects. 



5. Example: in case of k = 4 



In the following, we omit k = 4 in the subscripts of the notations. 



5.1. Example: random 4-noncrossing partitions. Set a^' 7 " denote the number of yQ 3 -walks 
of length s starting at (2, 1, 0), ending at (i, j, r) £ Q3 and set fl' 3 ' r = ail+i- 



Lemma 5.1. (a) Suppose 1 < I < N — 1 and (i,j,r) e W3, then 
(5.1) 

' _ f3<i-r _ fr,j.i _ fi,r,j , fj-r.i . rr.ij 

St J( JI Ji "r JI ' J£ 

f»^7> 1 fi—l,3,r , fi,j-l,r , fi,j,r-l _ f j,i,r _ f j—l,i,r f j,i-l.r _ f j,i,r-l 
J I ' J I "T 3 1 "T" J I J I J I J I J I 

_fr,j,i _ f r—l,j,i _ f r,j-l,i _ rr,j,i-l _ A,r,j _ f i-l,r,j _ A,r-l,j _ ri,r,j-l 

Ji Ji Ji Ji Ji Ji Ji Ji 



j,r—l,i 



+ 1 



j,r,i — l 



r—l,i,j 



+ 1 



r,i-l,j 



for s = 2i + 1, 



for s = 2£ + 2. 
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where for a, b, c e Z, /?' ,c = /or a + o + c>£ + 3 and /or £ = 0, 



I /or i = 2,j = 1, r = 0; 
otherwise. 



(b) ff' q ' s satisfies the recursion 

fP,q,s _ fp,q+l,s , rp+l,g,s . /.p,g,s+l »p— l,g+l,a . «>+l,g— 1,8 . /■p-l,g,s+l 
J£ ~~ J>-1 ->£-l J^-l ' J£-l H-l 

i /P+l,«,8— 1 i fP,9-l.s + l , y?p,9+l,s-l , x f p,q , /.p-l,(},s . f p,g-l : s . rp,q,s-l 
T "I- "I- "I- T "I- j£„! i" J 

Proof. We first prove assertion (a). Indeed, by definition of fg J ' r , we have 

(5.2) = and = + ff-l^ + fdi^ + /^T\ 

whence eq. (|5.1[) follows from Proposition 12.21 for the case k = 4. I.e. 
(0.3) oj e J = a e — a-i ~ a i — a i + a e + a t ■ 

Next, by construction, whence (b) and the proof of the lemma is complete. □ 



Once u)p J ' for arbitrary (i,j,r) £ W3 can be calculated follow the routine given in Lemma l5.ll we 
arrive at the following algorithm: 

Corollary 5.2. A random A-noncrossing partition can be generated, in 0(iV 4 ) pre-processing time 
and 0(iV 4 ) space complexity, with uniform probability in linear time. Each A-noncrossing partition 
is generated with 0(N) space and time complexity, see Algorithm 4- 



5.2. Example: 2-regular, random 4-noncrossing partitions. 

Lemma 5.3. The number ofB^, -walks ending at v £ W3 of length 11, where 1 < £ < N is given 
by 



(5.4) = £(-l) h '' 

h 



1,1, T 



h "2(e-h)+r 



Furthermore, we obtain the recurrence of cr^e-i given by 

fr. k\ i,j,r;* _ i,j,r;* , i-l,j,r;* , i,j-l,r;* , i,j,r-l;* 

{O.OJ °2l-l ~ u 2l-2 ' u 2t-2 ^ u 2t-2 * 21-2 
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Algorithm 4 Uniform generation of 4-noncrossing partitions 



Tableaux (Initialize the sequence of shapes, {X z }iZo N ) 
X° «- 0, X 2N «- 0, i «- 1 
while i < 27V do 
if i is even then 



X[0]<- 


- V(A* 


t l ,2N - 




1)) 


x[i]<- 


- V(Ai 




-(i- 


f 1)) 


X[2]<- 


- V(A* 




-(<- 


f 1)) 


X[3]<- 


- V(Ai 




-(i- 


f 1)) 


end if 











if i is odd then 

X[0]<- V(X l +\2N -(i + l)) 

X[2]^V(A!+ n 1 2 ,27V-(i+l)) 
X[3KV(A!g 3 ,2JV-(i + l)) 
end if 

sum <- X[0]+X[l]+X[2]+X[3] 

Shape <— Random(sum) (Random generates the random shape A^^ or _|-|^ with probability 
X[j]/ sum or A^ 1 with probability X[0]/sum) 
i<-i + l 

Insert Shape into Tableau (the sequence of shapes), 
end while 

Map( Tableau) (maps Tableau into its corresponding 4-noncrossing partition) 

Corollary 5.4. A random 2-regular, 4-noncrossing partition can be generated, in 0(N 5 ) pre- 
processing time and 0(N 4 ) space complexity, with uniform probability in linear time. Each 2- 
regular, 4-noncrossing partition is generated with O(N) space and time complexity, see Algorithm 
5, 6. 
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Algorithm 5 Uniform generation of 2-regular, 4-noncrossing partitions. 
Tableaux (Initialize the sequence of shapes to be a list {X l }\Zo N ) 



2 


A <- 0, A 


m «- 0, x 2 "- 1 «- 


0, i <— 1 








3 


while i < 2N - 1 do 










4 


if i is even then 










5 


X[0]<- 


V*(A l +1 ,27V- 


-1)) 








6 


x[i]<- 


V*(A^ i5 2iV- (i 


fi)),x[% 


- V*(A^,2iV- 


(< + 


1)) 


7 


X[3]<- 


V*(AV + i,2JV-(i 


H)),X[f 


-V*(Ai+ D \,2iV- 


(< + 


1)) 


8 


X[5]<- 


V^A^iV-^ 


+■ 1)),X[6]. 


- V*(A 4 l + ^,2iV- 


(< + 


1)) 


9 


else 












10 


X[0]<- 


W*(A^\2iV- (« 


f 1)) 








11 


x[i]<- 


W*(A!+Q i5 2iV- (i 


+ 1)),X[2]< 


- W^A^^iV 


~(i 


+ 1)) 


12 


X[3]<- 


W*(A!+^ 2 ,2iV- (i 


+ 1)),X[4]< 


- W*(A'j f p 3 ,2JV 


-(< 


+ 1)) 


13 


X[5]<- 


W*(A*+£,2i\T-(t 


+ 1)),X[6]- 


- W*(A!^ 3 ,2iV 


-(< 


+ 1)) 


14 


FLAG(flagO, flagl, flag2, flag3) 








15 


sum «- 


- X[0]+X[l]+X[2] + 


X[3]+X[4]4 


-X[5]+X[6] 







end if 

Shape <— Random(sum) (Random generates the random shape )?+Q.( or _Q.') with probability 
X[2j-l]/swm (or X[2j]/swm) or A^ 1 with probability X[0]/ sum) 
flagO <- 1, flagl <- 1, flag2 <- 1, flag3 <- 1 
if i is even and Shapc=A l f ^ i then 

flagl «- 
end if 

if i is even and Shapc^A 1 ^ then 

flag2 <- 
end if 

if i is even and Shape=A^p then 

flag2 «- 
end if 

if i is even and Shape=A^ 3 then 

flag2 «- 
end if 

if i is even and Shape=A^Tg then 

flag3 «- 
end if 

if i is even and Shapc=A^ 3 then 

flag3 «- 
end if 

if i is even and Shape=A t ~ 1 then 

flagO <- 
end if 

Insert Shape into Tableaux (the sequence of shapes) 
i <— i+ 1 
end while 

Map( Tableaux) (maps Tableaux into its corresponding 2-regular, 4-noncrossing partition) 
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Algorithm 6 FLAG: Distinguish the last odd step 
l: FL,4G(flagO, flagl, flag2, flag3) 
2: if flagO =0 then 
3: X[2], X[4], X[6]«- 
4: end if 

5: if flagl =0 then 

6: X[0], X[l], X[2], X[3], X[5]<- 

7: end if 

8: if flag2 =0 then 

9: X[l], X[2], X[3], X[4], X[5], X[6]<t- 
10: end if 

11: if flag3 =0 then 

12: X[l], X[3], X[5]<- 

13: end if 
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